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We have investigated noise-induced transition of atoms between double or triple phase-space 
attractors that are produced in the parametrically driven magneto-optical trap. The transition 
rates between two or three dynamic attractors, measured for various modulation frequencies and 
amplitudes, are in good agreement with theoretical calculations and Monte-Carlo simulations based 
on the Langevin equations. Our experiment may be useful to study nonlinear dynamic problems 
involving transitions between states far from equilibrium. 
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There have been much studies on fluctuation-induced 
transition between states that are in equilibrium or 
far from equilibrium. For equilibrium systems, after 
Kramers' seminal work |l|, many theories have been sug- 
gested and tested in many different experiments: for 
example, optically trapped Brownian particles 0-G1- ana- 
logue circuits p| , and semiconductor lasers @ . Recently, 
there have been several experiments on far from equi- 
librium systems, such as Penning trap Q, vibration- 
fiuidized granular matter Q, and Josephson junction 
Theoretical works on calculation of the transition rates 
[Tof or the transition paths , and analysis of the crit- 
ical exponents 0] have been performed. In particu- 
lar, interesting phenomena including oscillatory behavior 
with respect to the noise intensity [l|] or the phase [lil ]. 
and saddle-point avoidance have been also expected. 
However, except the Penning trap, the only quantita- 
tively investigated experimental system was the analogue 
electrical circuit |l6j| . which can be considered as ana- 
logue simulations. For the Penning trap, on the other 
hand, experiment was carried out within a very narrow 
parameter region: only near the bifurcation points. 

In this Letter, we report on the experimental study of 
the noise-induced transition in the paramterically driven 
magneto-optical trap (MOT) which is a far from equi- 
librium system. In particular, we have investigated the 
transition in nearly full parameter regions from super- 
critical bifurcation points (dynamic double attractors) 
to sub-critical bifurcation points (dynamic triple attrac- 
tors). For double attractors, we have measured the tran- 
sition rates by observing directly the change of popu- 
lation difference between two oscillating wells at vari- 
ous modulation frequencies and modulation amplitudes. 
On the other hand, the transition between three attrac- 
tors, which has never been probed experimentally to our 
knowledge, has been studied indirectly by measuring the 
populations in each well [ljj. We have also investigated 
the transition during one period of modulation and found 
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no phase dependence of the transition rates All 
the experimental results are found in good agreement 
with theoretical calculations and Monte-Carlo simula- 
tions based on the Langevin equations, which describe 
well the quantum dissipation 17] including spontaneous 
emission that is the origin of fluctuations of atomic mo- 
tion in MOT. 

As reported in previous papers |18| , when the intensity- 
modulation frequency / of the cooling laser is around 
twice the trap frequency between fx and ji, atoms are 
separated into two clouds [Fig. 1(a) Left] and oscillate in 
out-of-phase motion, which corresponds to the limit-cycle 
(LC) characterized by dynamic double attractors. Due 
to the nonlinearity of MOT, as / is increased above f% 
but below /a, there appears an additional stable attrac- 
tor at the center of the LC motion, called sub-critical 
bifurcation [Fig. 1(a) Right]. Given the modulation 
amplitude h, the corr esponding frequenci es are, fi = 
2 f - (/o/W/i 2 - h = 2/ o + (f /2)y/h?-h*, and 
fs = fo + (^/o/2/it)v4 + hj,, where /o is the natural 
trap frequency of MOT, and hx = P/i^fo) is the thresh- 
old value of h for parametric resonance to occur for a 
given damping coefficient f3. 

Assuming no fluctuations of atomic motion, initial con- 
ditions of atomic position and velocity determine which 
attractors (represented by red dots) an atom ends up 
with. Figure 1(b) shows each region of two attractors 
(Left) and that of three attractors (Right). For example, 
if the initial condition of an atom lies in the light gray 
region, the atom approaches the attractor located in the 
same basin. In reality, however, there exist atomic fluctu- 
ating motions due to spontaneous emission, resulting in 
broadened distributions of atomic position and velocity 
near the stable attractors. For large diffusion (or sponta- 
neous emissions), certain atoms may jump far from the 
original attractor and be transferred to another attrac- 
tor through the unstable regions near the boundary. In 
this case, the shape of the atomic phase-space distribu- 
tion resembles a dumbbell which has a narrow neck at 
the unstable point. The distribution represents the dy- 
namic potential wells, as discussed in the static potential 
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experiment 0]. 

The transition rate W n can be characterized by the ac- 
tivation law, W„ oc exp(— S n /D), which is the dynamic 
version of the Kramers' equation. Here S n is the acti- 
vation energy of each state n and D is the amount of 
diffusion which is proportional to the diffusion constant 
of MOT. In what follows we set n = 1 and 2 for the sta- 
ble states of two dynamic attractors and n — for the 
stationary state at the MOT center. As observed in the 
stationary position-space potentials, the transition rate 
increases as the diffusion increases or the activation en- 
ergy decreases. Note that since all the nonlinear terms in 
the equations of motion of MOT have the same sign as 
the first order term, the activation energy is proportional 
to / and h. 

Unlike the single electron in the Penning trap, there 
are more than 10 7 atoms in the initial MOT and each 
cloud of the LC motion under parametric excitation is 
almost equally populated. In order to monitor the tran- 
sition between two dynamic attractors, it is needed to 
blow away one of the clouds because the atomic number 
transferred from attractor 1 to 2 is the same as that from 
2 to 1. When atoms, say in attractor 1, are removed, one 
can observe some atoms in attractor 2 are transferred to 
1 with time. Due to the two-way transitions the popula- 
tion difference decreases exponentially. For many particle 
systems the transitions can be described by the following 
simple rate equations, 



dNi 

~dT 

dN 2 

dt 



R-ryNt- W 1 N 1 + W 2 N 2 , 



R-jN 2 - W 2 N 2 + WxN! , 



(1) 



where Ni (N 2 ) is the population of the attractor 1 (2), J? 
is the trap loading rate, 7 is the loss rate due to collisions 
by the background atoms, and W\ (W 2 ) is the transition 
rate from the attractor 1 to 2 (2 to 1). Since R and 
7 are small in our experiment (as discussed later), we 
neglect these terms and also assume W — W-y — W 2 due 
to the symmetry. Then the steady-state solution of the 
population difference AiV(= N 2 — Ni) is given by 
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where AN® is the difference with one attractor (say, n = 
1 state) empty, or N% = 0. 

Our experimental setup is similar to that described 
in previous work |l8j], but with new features: we have 
used a photodiode array 0] to measure the transition 
rates and used a resonant laser to blow away selectively 
one atomic cloud (say, n = 1). The blowing laser was 
cylindrically focused at 5 mm left from the center of the 
LC motion. The intensity of the laser was over 5 times 
saturation intensity and it was turned on for 3 ms in 
order to remove only one atomic cloud. 

Experimental results for double (triple) attractors are 
obtained at the intensity of 0.034 I s (0.06 I s ) and the 
detuning of -2.7 T for the cooling laser along the longitu- 
dinal z-axis, and the magnetic field gradient of 10 G/cm, 
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FIG. 1: (a) The atomic fluorescent images of the LC mo- 
tion when there are two attractors (Left) and three attractors 
(Right), (b) The phase-space map of two (Left) and three 
attractor (Right) conditions. 
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FIG. 2: (a) Typical contour plot of absorption signals ob- 
tained at / = 63 Hz (f/f = 1.9) and h = 0.7. (b) Decay of 
the population of state 2 (filled box), growth of the atomic 
number of the empty state 1 (empty box) as measured from 
(a), and the averaged half of the total number (asterisk), (c) 
The population difference between the two states shows ex- 
ponential decay. 



where I s = 3.78 mW/cm 2 is the averaged saturation in- 
tensity and r = 2-7TX 5.9 MHz is the natural line width. 
The laser intensity and detuning along the transverse 
x- and y-axis are 0.63 I s and -3.0 T, respectively. For 
these parameters, the measured trap frequency is 33.4 
Hz (43.9 Hz, for triple states), and the detected damping 
coefficient is 45.4 s _1 (for triple states, 85.9 s _1 ), which 
are in good agreement with the simple Doppler theory. 
Note that the sub-Doppler nature of MOT is dramati- 
cally suppressed and thus neglected when the transverse 
laser detuning is slightly different from that of the z-axis 
laser, as discussed in depth in our previous work [Tsj . 
The typical data about the transition between dynamic 
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double wells are presented in Fig. 2. Figure 2(a) shows 
a contour plot generated by the atomic absorption sig- 
nals recorded on the 16 channel photo-diode array. The 
vertical axis represents the atomic position and the lon- 
gitudinal axis the time evolution. In Fig. 2(a), large ab- 
sorption (or large atomic number) is indicated by bright 
color. From the plot one can trace the oscillating LC 
attractors and measure the atomic number. It can be 
observed that one attractor is made empty at 7 half cy- 
cles from the left, which is then rcpopulatcd with time, 
as plotted in Fig. 2(b). Figure 2(c) presents the tem- 
poral variation of the population difference (AN), which 
clearly shows the exponential decay as expected in Eq. 
J2J where the decay rate is twice the transition rate. Note 
that the total number of atoms is nearly conserved dur- 
ing the transitions [asterisks in Fig. 2(b)], which justifies 
neglect of the loading rate R and loss rate 7 in Eq. 
In reality, the loss rate of our MOT is 0.15 s" 1 which is 
much smaller than the typical transition rates. 
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FIG. 3: Transition rates for various (a) modulation fre- 
quencies and (b) modulation amplitudes, compared with the 
Monte-Carlo simulations for scaled diffusion D. 

The experimental results in Fig. 3 were obtained by 
exponential fitting of Fig. 2(c) for (a) various modula- 
tion frequencies at ft, = 0.72 and (b) modulation ampli- 
tudes at f/fo — 1-89. The curves in Fig. 3 represent 
the Monte-Carlo simulation results obtained from simple 
Doppler equations based on two-level atom and random 
spontaneous emission. We have first calculated the tra- 
jectories of atomic motion of 10 4 atoms initially in n = 2 
state and monitored the subsequent time evolution of the 
population of the two states. We have then obtained the 
transition rates by fitting Eq. @ ■ As can be observed in 
Fig. 3, the experimental results are in good agreement 
with the simulations when diffusion constant is 0.7 times 
DD, where DD represents the Doppler diffusion coeffi- 
cient (Eq. (22) in Ref. |U). 

In the low damping regime, one can also find that the 
transition rates approximately vary around the bifurca- 
tion point as 



W = c\ exp 



hD U 



(3) 



where /i is the modulation frequency at the bifurcation 
point, and ci and C2 are constants |10j . The constant 
C2 = (l//o) y/Z/hr was obtained by direct calculation 



[To| . The constant c\ and the diffusion D were used as fit- 
ting parameters. The consequent diffusion D for our ex- 
perimental conditions is 0.071(8), which are in the same 
order to the calculated value with simple Doppler theory 
D = 3A / (2m 2 uj f3 2 ) DD [hJ. Here the coefficient A a is 
nonlinear coefficient of a term z 3 in the equation of mo- 
tion of a parametrically-driven atom, which is presented 
as Eq. (7) in Ref. [gof. The triangles in Fig. 3(a) and 
3(b) represent the fitted results. Even though the calcu- 
lations are mainly well applied in the limit / j> /j and 
/3 < wo, we can see these fittings very well describe the 
experimental results. 

Since W and / are of the same order in our experi- 
ment, we also have studied the time dependence of the 
transition during one period of parametric oscillation. 
By averaging the absorption signals over 20 times, we 
have obtained the population evolution during each pe- 
riod. The experiments were performed at / = 1.76/o an d 
h = 0.72, with the same laser conditions and magnetic 
field gradient as above. The transition rate W was 15.4 
s , which was measured during three oscillation periods 
until the population difference became negligible. In par- 
ticular, we have not found any phase dependence of the 
transition rates, which indicates that transitions occur 
constantly without any dependence of location of atomic 
clouds in phase space. Note that these results are very 
different from those observed in periodic driven systems 
[Til , which were studied only theoretically. Experimental 
results are also in very good agreement with simulations, 
which do not show any clear variation of the transition 
rates during one period. 
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FIG. 4: (a) Contour plot of the absorption signals for triple 
well conditions, where /=106 Hz (= 2.42 / ) and h=0.72. (b) 
Contour plot of the simulation results, where /=105.3 Hz (= 
2.41 /o) and /i=0.7. (c) Relative populations of state 1 and 2 
(iVi + N 2 )/N T compared to that of state (N /N T ), where 
points (filled box and filled circle) represent the experimen- 
tal results and curves the simulations. Filled triangles show 
the experimental ratios of Ni + N2 to No and curves the 
simulation results. 

For the case of dynamic triple attractors, direct obser- 
vation of the transition is difficult in our setup. When 
!i < / < /3j h is numerically observed that transition 
occurs only between one of the two dynamic states (1 
or 2) and the stable state (0), with no direct transitions 
from 1 (2) to 2 (1). Moreover, unlike the double wells, 
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the transition rates Wio from state i to are not nec- 
essarily the same as Woi from state to i (i = 1,2). 
Thus the method of just removing the state 1 and mon- 
itoring the time evolution of the population difference 
does not provide the necessary information of the rel- 
evant transition rates. Therefore we have studied the 
transition in the sub-critical region by using an indirect 
method of measuring the populations in each state. The 
ratios of the population Ni (or N 2 ) to Nq can be used 
to obtain the transition rates, that is, Ni/N = N 2 /N 
=Wq/2W, where we have employed W w — W 20 = W 
and W01 — Wq 2 — Wo/2. Note that for a larger modula- 
tion frequency, N becomes larger than Aq (or N 2 ) |1C| . 
which is confirmed in our experiment. 

The experimental and simulation results for dynamic 
triple attractors are presented in Fig. 4. Figure 4(a) 
shows the motions of atomic clouds during single pe- 
riod of oscillation, which provides the population in each 
state. States 1 and 2 show out-of-phase oscillating mo- 
tions, while state stays at the center. Our simulations 
[Fig. 4(b)] have very similar results to the absorption im- 
age [Fig. 4(a)], with different spatial resolutions (1 mm 
for Fig. 4(a) and 0.1 mm for Fig. 4(b)). 

The data of (N t + N 2 )/N T (filled boxes) and N /N T 
(filled circles) in Fig. 4(c) are obtained by fitting the 
profiles of Fig. 4(a) at 0.58 it phase. The dashed curves 
{{Ni + N 2 )/N T ) and dotted curves (N /N T ) are derived 
from Fig. 4(b) at 0.5 n. Here N n (n = 0,1,2) is the pop- 
ulation of the state n and Nt—Nq+Ni+N 2 . The profile 
shows three peaks, representing three states, are fitted by 



triple Gaussian functions. Filled triangles (experiment) 
and solid curves (simulation) are just the ratios of N\ + 
N 2 to Nq, which also provides the ratio of the transitions 
rates (Wq/W). One can observe No is increased as / 
increases, whereas N\ and N 2 are decreased. It shows 
that the transition from state 1 (2) to state becomes 
easier than that from to 1 (2) as / increases. These re- 
sults are in good agreement with theoretical calculation 
Wo/W = exp[(Si — So)/D] where S n is the activation 
energy of state n [To| . 

The parametrically driven MOT is an ideal experimen- 
tal realization of dynamic double and triple attractors 
under the conditions far from equilibrium. We have in- 
vestigated the noise-induced transition at various / and h 
from the super-critical to sub-critical bifurcation region. 
The results are well described by the Monte-Carlo simu- 
lations based on the Langevin equations and the simple 
Doppler theory. Note that our system may be important 
and useful to study nonlinear problems associated with 
transitions in far from equilibrium systems. One may 
also extend our system to the study of transitions be- 
tween dynamic states due to non-symmetric oscillations 
that occur in more nonlinear situations |2lj . 
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